Draft version February 2, 2008 

Preprint typeset using I^T^ style emulateapj v. 04/03/99 



PREDICTING PLANETS IN KNOWN EXTRA-SOLAR PLANETARY SYSTEMS I: TEST PARTICLE 

SIMULATIONS 

RoRY Barnes', Sean N. Raymond' 



E-mail:rory@astro. washington.edu 

Draft version February 2, 2008 



o 
o 

(N 

X) 

(D 

vn 

(N 

> 
(N 

(N 

o 
o 

Oh 

I 

o 



ABSTRACT 

Recent work has suggested that many planetary systems lie near instability. If all systems are near instability, an 
additional planet must exist in stable regions of well-separated extra-solar planetary systems to push these systems 
to the edge of stability. We examine the known systems by placing massless test particles in between the planets 
and integrating for 1-10 million years. We find that some systems, HD168443 and HD74156, eject nearly all test 
particles within 2 million years. However we find that HD37124, HD38529, and 55Cnc have large contiguous 
regions in which particles survive for 10 million years. These three systems, therefore, seem the most likely 
candidates for additional companions. Furthermore HD74156 and HD168443 must be complete and therefore 
radial velocity surveys should only focus on detecting more distant companions. We also find that several systems 
show stable regions that only exist at nonzero eccentricities. 

Subject headings: 



I. INTRODUCTION 

Several categories of planetary systems have been discovered 
in the past several years. Some appear dynamically similar to 
our own Solar System (SS), but others appear quite different 
(Barnes & Quinn 2004, hereafter BQ). Whether these systems 
really fall into unique categories, or are the first examples of a 
continuous spectrum of stability remains to be seen. Nonethe- 
less we will define three categories of systems: resonant, in- 
teracting, and separated. For the purposes of this paper, we 
will define a planetary system as a star with 2 or more com- 
panions; at least a 3-body system. Resonant systems contain 
2 or more planets in mean motion orbits. The GJ876 (Butler 
et al. 2001) and HD829432 systems are in 2:1 resonance, and 
the 2 inner planets of the 55Cnc system (Marcy et al. 2002) 
are in 3:1 resonance. The interacting systems of v And (But- 
ler et al. 1999), HD 12661 (Fischer et al. 2003), 47UMa (Butler 
et al. 2001), and the SS are not in resonance, but their orbits 
are close enough that the planets may perturb each other. The 
final category is the separated systems. In these systems, the 
planets are separated enough that they are not interacting on 
long timescales (^ 10^ years). This paper will focus on these 
separated systems, specifically HD 168443 (Marcy et al. 2001), 
HD74156^ HD38529 (Fischer et al. 2003), and HD37 124 (But- 
ler et al. 2003). We will also examine 55Cnc, which is a com- 
bination of a resonant and a separated system. 

Studies of the known resonant and interacting systems have 
demonstrated that these systems are on the verge of instability 
(BQ). Slight changes in orbital elements, specifically eccentric- 
ity or proximity to perfect resonance, lead to a catastrophic dis- 
ruption of the systems. Planetary systems near instability are 
as tightly packed as possible; there is no room for additional 
companions. In this paper we explore the possibility that all 
systems are tightly packed, implying that additional compan- 
ions lie between the extant planets in separated systems. In this 
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"packed planetary systems" (PPS) hypothesis, the undetected 
planets have not been observed because there are not enough 
data to discover the additional planet, or the planetary mass 
falls at or below the detection limit of current Doppler technol- 
ogy. BQ also examined the stability of the Sun-Jupiter-Saturn 
system and showed that it lies further from instability than the 
complete gas giant system. The Sun-Jupiter-Saturn model is 
the only system known to be incomplete and therefore supports 
the hypothesis that all systems lie near instability, and are hence 
packed. 

The PPS model assumes that planet formation is an efficient 
process. As many planets (or at least gas giants) form in a cir- 
cumstellar disks as possible. This is seen in the gas giant region 
of the SS, as well as in the resonant and coupled systems. In this 
scenario, as dust congeals into ever larger bodies, the mutual 
gravitational forces between protoplanets perturb each other in 
ever increasing amounts as the masses grow. As the protoplan- 
ets become larger, they acquire the ability to pull other planets 
into favorable positions (dynamically stable), or, if the pertur- 
bations become strong enough, eject them completely from the 
system. The implication is that any stable region in between 
known planets harbors an (as yet) unseen companion. In this 
paper we attempt to map out stable regions in semi-major axis 
and eccentricity space of known separated systems. 

We perform numerical simulations using test (massless) par- 
ticles to search for stable regions in HD168443, HD74156, 
HD37124, HD38529, and 55Cnc. In addition to these systems 
we also perform a control experiment. We replace the middle 
planet of the v And system with a belt of test particles to evalu- 
ate how well this experiment predicts the orbit v And c. We find 
that some systems, HD74156 and HD168443, contain no stable 
regions and are complete to the most distant planet. HD37124, 
HD38529 and the resonant/separated system 55Cnc contain a 
zone of stability which may harbor an undetected planet. This 
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paper is divided into the following sections. In §2 we describe 
the numerical techniques used to examine these systems. In §3 
we present the results for systems which may have additional 
companions. In §4 we present the results of the v And system, 
with the middle planet removed. We summarize our results in 
§5 and suggest future work in this field. This is the first in a 
series of papers which will explore the possibility of additional 
planets in known systems. In Paper II we will place massive 
planets in the stable regions identified by this work to refine 
our predictions. In Paper III we examine how planet formation 
might occur between the detected planets. 

2. NUMERICAL METHODS 

These simulations were performed with SWIFT"* (Levison & 
Duncan 1994). Specifically we used the regularized mixed vari- 
able symplectic integrator, RMVS3. This code is designed to 
quickly integrate a system through close approaches (via regu- 
larization), yet still maintain conservation of energy (via a sym- 
plectic algorithm). For most systems (except 55Cnc and v And) 
we integrate for 10^ years. We require every simulation to con- 
serve energy to a factor of lO""*, which has been shown to be ac- 
curate enough for simulations of planetary system stability on 
these timescales (BQ). Generally we conserve energy to better 
than 1 part in 10^. 

The systems presented here are all coplanar. The mass- 
inclination degeneracy of Doppler observations is broken by 
assigning the companions' masses to be the observed mini- 
mum masses. All other orbital elements are their best fit val- 
ues. Note that in coplanar systems, the longitude of ascending 
node is meaningless. Test particles are all also coplanar. They 
are spaced every 0.002AU in semi-major axis space, and every 
0.05 in eccentricity space. Their mean anomalies are placed 
uniformly throughout [0,27r). In non-circular trials the lon- 
gitude of periastton of the test particles was aligned with the 
most massive planet. This alignment is arbittary, but probably 
valid for low e (< 0.1). However some research suggests that 
for e > 0.1 the apses should have been anti-aUgned (Laughlin, 
Chambers & Fischer 2002), although this phenomenon is not 
strictiy necessary (BQ). 

This procedure is similar to other work on test particles in 
V And (Rivera & Lissauer 2000), GJ876 (Rivera & Lissauer 
2001), and 55Cnc (Rivera & Haghighipour 2002). The major 
difference being that we examine non-circular test particle or- 
bits. Additionally Menou & Tabachnik (2003) placed test parti- 
cles in the habitable zones of these systems. With the exception 
of HD37124, we recover their results. 

3. RESULTS 

In Table 1, the orbital elements of all the planetary systems 
used in this paper are shown. As there are no published data 
for some systems, we are forced to use data from the inter- 
net, which are constantly changing. In §3.1 we examine the 
complete systems of HD74156 and HD 168443 and in §3.2 we 
analyze the HD37124, HD38529, and 55Cnc systems which 
contain zones of stability. 

3.1. Complete Systems: HD74156 & HD168443 

From Table 1 nothing about these 2 systems would indicate 
a priori that they would have no stable zones, although their 
planetary masses are some of the largest of this subset of sys- 
tems. Note, however, that HD38529c is the 2nd largest planet 
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examined, and it has one of the largest zones of stability (see 
§3.2). The ratios of the periods, R, of HD74156 and HD168443 
are 51.4 and 30.5, respectively, which are not the smallest val- 
ues among separated systems. Nonetheless these systems show 
littie evidence for stable particles between the currentiy known 
planets. 

In HD168443, /?=30.5, and we need 786 test particles to fill 
the region between planets b and c. In Figure 1, we see that 
regardless of eccentricity or semi-major axis no test particles 
survive for even 2 million years. We therefore conclude that 
this system is complete out to planet c, there can be no aster- 
oid belt in this system, and radial/astrometric surveys should 
sample this system sparsely. 




Fig. 1.— Stability of test particles in HD168443. The height of the 
bars corresponds to how long test particles remained bound to the par- 
ent star. Regardless of eccentricity, no test particles can survive in this 
system for even 2 milUon years. 

The best fit orbital elements for HD74156 have fluctuated 
throughout the past several years. The elements presented in 
Table 1 are from the discovery website and date to May of 2002. 
The elements changed dramatically in August of 2002 (namely 
flc increased to 3.8AU). The data have recently been published 
(Naef et al. 2003), and the elements changed to values similar 
to those of May 2002. Although the system presented here is 
slightly different than the current best fit, our work is similar 
enough that any differences should be negligible. 

HD74156 consists of two planets separated by nearly a factor 
of 40 in period. We need 822 test particles for this system. The 
eccentricities and masses are lower than in HD 168443. In Fig- 
ure 2, we present the survivability of test particles in this sys- 
tem. A very narrow band of stability may exist at 0.6AU, this 
strip is likely unstable on timescales larger than 10^ years. Note 
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Table 1 - Orbital Parameters of Selected Planetary Systems 



System 


Planet 


M (Mj) 


a (AU) 


e 


VD 


T (JD) 


HD168443 


b 


7.73 


0.295 


0.53 


172.9 


2450047.58 




c 


17.2 


2.87 


0.2 


63 


2450250.6 


HD74156 


b 


1.56 


0.28 


0.649 


183.7 


2451981.40 




c 


7.5 


3.47 


0.395 


240.0 


2450849.0 


HD37124 


b 


0.86 


0.1 


0.54 


97.0 


2451227 




c 


1.01 


0.4 


2.95 


265.0 


2451828 


HD38529 


b 


0.78 


0.129 


0.29 


87.7 


2450005.8 




c 


12.8 


3.68 


0.36 


14.7 


2450073.8 


55Cnc 


b 


0.84 


0.115 


0.02 


99.0 


2450001.479 




c 


0.21 


0.241 


0.339 


61.0 


2450031.4 




d 


4.05 


5.9 


0.16 


201.0 


2452785 


V And 


b 


0.69 


0.059 


0.015 


32.0 


2450001.588 






1.96 


0.83 


0.25 


251.0 


2450160.1 




d 


3.98 


2.50 


0.34 


255.0 


2450044.0 



^Not used in these simulations 
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that this strip only develops at e > 0.1. Therefore any asteroids 
that may exist in this system must be on significantly eccentric 
orbits. Although it appears unlikely that any planet could exist 
in the putative region, we will show in Paper II (Raymond & 
Barnes, in preparation) that Saturn mass companions can sur- 
vive in the system for at least lO^yr. 




Fig. 2. — Stability of test particles in HD74156. The stable region 
located at a = 0.6Ai7,0.1 < e < 0.25 is so narrow that an additional 
companion is very unlikely. Note that the band arises at nonzero ec- 
centricity. 



3.2. Candidate Systems: HD37124, HD38529, & 55Cnc 

Three systems show broad regions of stability for test parti- 
cles. The current orbital parameters of HD37124, HD38529, 
and 55Cnc are presented in Table 1. HD37124 and HD38529 
are classical separated system, whereas 55Cnc is a resonant sys- 
tem with a distant, separated companion. In 55Cnc the outer 
two planets are separated by R=\2\, HD37124, /?=10.1, and 
HD38529,/?=152. 

In Figure 3, we show the stable zone for HD37124. For this 
system R=\2.1 and we integrated 588 test particles. Although 
dotted with ejections due to high order resonances, there do 
appear to be significant regions of stability, most notably is at 
a semi-major axis of approximately lAU. There is also an in- 
crease in stability near 0.15 < e < 0.2. As in previous systems, 
we note that the stable region is not largest at e=0. Therefore 
we suggest that the most likely orbit for an additional compan- 
ion is at a semi-major axis slightly smaller or larger than lAU 
(which corresponds to the 5:2 resonance with planet b), and an 
eccentricity near 0.15. Should there be no planet here, these re- 
sults suggest that the presence of an asteroid belt in this region 
of phase space is likely. 




Fig. 3.— Stability of test particles inHD37124. The most likely orbit 
for an additional companion lies just interior or exterior to lAU, with 
an eccentricity of ~ 0.15. This plots also shows the possible orbits of 
asteroids should no additional planet be present. 

Note also in Fig. 3 that an additional band of stability arises 
at a « 1 .2AU above e=0. 1 . A third band is visible at a = 0.9AU, 
e = 0.15. A fourth band arises at a « 1.7AU at e > 0.3. This 
suggests that the significant {i.e. larger than 0.1) values of extra- 
solar planetary eccentricities may not be detrimental to system 
stability. This has been suggested before (Murray & Dermott 
1999, Menou & Tabachnik 20003), but is shown dramatically 
in Fig. 3. Although moderate eccentricities may encourage sta- 
bility, at some point in every system, a critical eccentricity is 
reached and the system becomes unstable. For HD37124 this 
threshold is at e=0.3. 

As mentioned in §2, Menou & Tabachnik obtained a differ- 
ent result for this system. They integrated the region from 0.6 to 
1.2AU for 10* years, with test particles separated by 0.006AU. 
They found that no test particles, with eccentricities close to 
0.05, survived in this region. The discrepancy results from 
their stringent criteria for ejection, such as labeling test parti- 
cles which cross the boundary of the habitable zone as unsta- 
ble. We find test particles in our habitable zone (in semi-major 
axis space) often have significant eccentricities, and most likely 
crossed the boundary at some point during the simulation. 

In Figure 4 the stable regions of HD38529 are shown. R 
for this system is 152, and hence we integrated 1092 test par- 
ticles. This region also contains a narrow zone of stability. 
In this system we see that stability lies in between 0.25 and 
0.75AU, and e < 0.3. The interior edge of this stable zone is 
quite sharp, probably because the mass of the inner planet is 15 
times smaller than the outer. This system has a larger zone of 
stabiUty than HD37124, so it seems more Ukely that a planet 
nught exist in this system. 
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Fig. 4.— Stability of test particles in HD38529. Stability is most 
likely in the region 0.25AU <,a <^ 0J5AU. Above e = 0.3 the proba- 
bility of survival decreases dramatically. 



Fig. 5. — Stability of test particles in 55Cnc. A 1.9AU wide region 
of stability exists in this system. Unlike other systems, high eccentric- 
ity does not appear to promote stabiUty. This system seems the best 
candidate for additional companions. 



55Cnc is the only system with a resonant pair, and a distant 
companion. Of all the systems examined here, R is the largest. 
Because of this wide separation, 55Cnc requires 2164 test par- 
ticles and we could only integrate the system for 5 milhon years 
due to resource constraints. In Figure 5 we plot survivability in 
this system. Not surprisingly, this system shows the broadest 
range of stability for test particles. For all eccentricities there 
appears to be stabiHty from 0.9AU to 2.8AU. From 2.8AU to 
3.5AU the stable regions are broken up by mean motion reso- 
nances with planet d. There does appear to be a sUght prefer- 
ence for low eccentricities in this system. Therefore we suggest 
that a planet might exist close to 2AU (the center of the stabil- 
ity zone), and with e ^ 0. 1 . As before we also suggest that this 
region might harbor an asteroid belt or terrestrial planets if it 
does not contain a gas giant. 

55Cnc was also examined by Rivera & Haghighipour (2002). 
They examined test particles between planets c and d for 5 mil- 
lion years, but all the test particles had zero eccentricity. They 
find a similar stable region. This is also consistent with the 
results of Marcy et al. (2002) which stated that an additional 
planet at lAU would be stable. Rivera & Haghighipour also 
examined the region exterior to planet d and found stability did 
not appear again until beyond lOAU. Given the similarity be- 
tween our results and theirs, it is Ukely that we would obtain a 
similar result, even for nonzero eccentricities. 



Although there is a broad region of stability, this is only true 
of massless particles. The resonant pair in this system hes close 
to instabiUty (Henderson & Bames, in preparation), therefore 
an additional perturbing mass may increase the eccentricity, or 
move the system from perfect resonance, and hence destabilize 
the system. In Paper 11 we will address this issue when we place 
massive companions in this region. 



4. A CONTROL: v ANDROMEDAE 

To test our methodology we simulated v And to determine 
if we could predict where the middle planet (planet c) might 
he. In Table 1 we presented the orbital elements of the planets 
in this system. In our 3 body model, /?=284 and we required 
814 test particles. The region between planets b and d proved 
so stable that we were only able to integrate this system for 1 
million years due to CPU time limitations. The stability of this 
system is shown in Figure 6. As we can see from this figure, 
we do indeed find that planet c hes in the stable region. 
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Fig. 6. — Stability of test particles in v And. This system shows 
stability between 0.1 and 1.2AU. Note also that the system is more 
stable at eccentricities larger than 0.2. The current values of Uc and Cc 
are 0.8AU and 0.24, respectively. 

Test particle simulations of v And were also performed by 
Rivera & Lissauer (2000). Their simulations ran for 5 million 
years. They used all three planets, but did place test particles 
between planets b and c, in addition to particle beyond the orbit 
of planet d. They found that a very narrow region of stability be- 
tween planets b and c of approximately 0.35AU in width. They, 
too, find a sharp edge near planet b. Our hypothesis should 
therefore suggest that an additional planet may lie in this re- 
gion between planets b and c as well (see §5). However this is 
the narrowest zone of stability seen in any system, and, given 
that V And already lies near the edge of stability, an additional 
planet located here may disrupt the system. 

5. DISCUSSION AND CONCLUSIONS 

Building on the work of BQ we have argued that additional 
planets should exist in separated systems so that they, too, lie 
close to instability. We have tested the PPS hypothesis by inte- 
grating the orbits of a large number of massless test particles in 
five known extra-solar planetary systems. The results of these 
simulations are summarized in Table 2. In this table A'^ is the 
number of test particles used, Aa and Ae are the approximate 
zones of stability where an additional companion might exist. 
We find that some systems cannot contain additional planets 
(HD74156 and HD168443), while others have significant sta- 
ble zones (HD37124, HD38529, and 55Cnc). From our control 
experiment, we see that Figure6 most resembles FigureS, again 
suggesting that 55Cnc is the most likely candidate for an addi- 
tional companion. 

We must note one ambiguity in the PPS hypothesis. As has 
been shown in other work (Barnes & Quinn 2001, BQ), v And 



is already near the edge of stability. Perhaps an additional com- 
panion was in this system, but it was dynamically unstable, and 
ejected. This left c and d on interacting orbits, and a small an- 
nulus of stability. Planetary systems with any 2 planets on the 
edge of stability may negate the possibility of predicting addi- 
tional companions. They are already on the edge, and hence 
need no additional companions to push them there. This is es- 
pecially relevant in 55Cnc. The resonant pair is on the edge 
(Henderson & Barnes, in preparation), so perhaps the large gap 
between planets c and d is irrelevant. However there is an ad- 
ditional complication in comparing v And and 55Cnc. Namely 
they may have completely different formation histories. Reso- 
nant systems most likely formed from resonant capture during 
the migration epoch (Snellgrove, Papaloizou, & Nelson 2001), 
whereas v And may have formed from a large scattering event 
(Rasio & Ford 1996, Malhotra 2002, BQ). If a planetary system 
needs only 2 planets near the edge to be "packed", then 55Cnc 
is already packed, and has no additional companions. 

Here we have suggested that additional planets exist in be- 
tween the known planets. It could be, though, that additional 
planets Ue beyond these in an interacting configuration. The 
work of Rivera & Lissuaer (2000) shows that this is the case for 
the V And system (see §4). If the separated systems presented 
here actually have more distant companions on interacting or- 
bits, these systems would still lie close to the edge. Examina- 
tions of this possibiUty, however, were beyond the scope of this 
paper. 

We have also found that some test particles are more stable 
on significantly eccentric orbits. This is true for HD74156 and 
HD37124. Although there has been some mention of this phe- 
nomenon {i.e. Menou & Tabachnik 2003), a satisfactory expla- 
nation for this observation is not apparent. Future work should 
address this issue. 

To further explore the PPS scenario, we will integrate the full 
4 or 5 body systems with a massive planet in the regions of sta- 
biUty defined here. This work will be Paper II. Additionally 
we will explore the formation of terrestrial planets in these sys- 
tems in Paper III. In these papers we will show more evidence 
that the predicted planets can survive for upward of 100 milhon 
years. 

Of course stable regions do not guarantee stability, but it 
is nonetheless intriguing to postulate their existence. If the 
candidate systems do contain additional companions, then we 
strengthen the theory that planet formation is an efficient pro- 
cess and add this as a new requirement to models of planet for- 
mation. Whether or not this may break the degeneracy between 
the so-called core accretion model (i.e. Pollack et al. 1996; Bo- 
denheimer, Hubickyj, & Lissauer 2000) and the gravitational 
coUapse model {i.e. Boss 2002; Mayer et al. 2002) remains to 
be seen. 

Should the planets predicted here be discovered, it would 
mark the first prediction of a planet since Neptune by John 
Couch Adams in the 19th century. The radial velocity surveys 
are providing a rich source of knowledge for the field of planet 
formation, altering our understanding of our own solar system, 
and changing perspectives of our place in the universe. The 
correct prediction of a new planet would represent a major step 
toward understanding the mechanisms of planet formation. 
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Table 4 - Results of Test Particle Simulations 



System 


R 


N 


Aa 


Ae 


HD168443 


30.5 


786 






HD74156 


51.43 


1004 






HD38529 


152.0 


1092 


(0.27,0.82) 


(0.0,0.3) 


HD37124 


12.7 


588 


(0.895,1.11) 


(0.0,0.25) 


55Cnc 


121.0 


2164 


(0.72,2.77) 


(0.0,0.25) 


V And 


284 


814 


(0.1,1.2) 


(0.0,0.3) 
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